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Abstract 

A new class of linear second order hyperbolic partial differential operators satisfying 
Huygens' principle in Minkowski spaces is presented. The construction reveals a direct 
connection between Huygens' principle and the theory of solitary wave solutions of the 
Korteweg-de Vries equation. 



Mathematics Subject Classification: 35Q51, 35Q53, 35L05, 35L15, 35Q05. 



1 



I. Introduction 



The present paper deals with the problem of describing all linear second order partial differential 
operators for which Huygens' principle is valid in the sense of "Hadamard's minor premise". 



Originally posed by J.Hadamard in his Yale lectures on hyperbolic equations this problem is 
still far from being completely solvedQ 

The simplest examples of Huygens' operators are the ordinary wave operators 

in an odd number n > 3 of space dimensions and those ones reduced to (|l]) by means of elementary 
transformations, i.e. by local nondegenerate changes of coordinates x ^-^ f{x) ; gauge and 
conformal transformations of a given operator C ^ 9{x) o £, o 6{x)~^ , C ^ n{x)C with some 
locally smooth nonzero functions 9{x) and fi{x). These operators are usually called trivial Huygens' 
operators, and the famous "Hadamard's conjecture" claims that all Huygens' operators are trivial. 

Such a strong assertion turns out to be valid only for (real) Huygens' operators with a constant 
principal symbol in n = 3 |Q . Stellmacher |^0| found the first non-trivial examples of hyperbolic 



wave-type operators satisfying Huygens' principle, and thereby disproved Hadamard's conjecture 
in higher dimensional Minkowski spaces. Later Lagnese & Stellmacher |^ extended these ex- 



amples and even solved Hadamard's problem for a restricted class of hyperbolic operators, 
namely 

£ = □„+! + , (2) 

where u is an analytic function (in its domain of definition) depending on a single variable 
only. It turns out that the potentials u{z) entering into (H) are rational functions which can be 
expressed explicitly in terms of some polynomials^ Vk (z) : 



u{z) = 2\^—j \ogVkiz) , A; = 0,l,2,..., (3) 

the latter being defined via the following differential-recurrence relation: 

rUiVk-i - n-i^+i = + l)Vl , Vo = l,V, = z. (4) 

Since the works of Moser et al. @], @ the potentials (^) are known as rational solutions of the 
Korteweg-de Vries equation decreasing at infinity^. 

"'^Hadamard's problem, or the problem of diffusion of waves, has received a good deal of attention and the 



literature is extensive (see, e.g., [[12|, |15[, ]21| 



35[ , and references therein) . For a historical 



account we refer the reader to the articles lit 

^This remarkable class of polynomials seems to have been found for the first time by Burchnall and Chaundy 

[0. 

The coincidence of such rational solutions of the KdV-hierarchy with the Lagnese-Stellmacher potentials has 
been observed by Schimming pm, [p9|. 
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A wide class of Huygens' operators in Minkowski spaces has been discovered recently by Veselov 
and one of the authors 0, |T^ (see also the review article [^). These operators can also be 



presented in a self-adjoint form 

£ = + uix) (5) 

with a locally analytic potential u [x) depending on several variables. More precisely, u (x) belongs 
to the class of so-called Calogero-Moser potentials associated with finite reflection groups [Coxeter 
groups) : 

Em„(ma + l)(a,a) 
• (6) 



U(X) , , _ 

[a, xy 



In formula (|) 3?+ = 3?+(^) stands for a properly chosen and oriented subset of normals to 
reflection hyperplanes of a Coxeter group Q. The group Q acts on M"+^ in such a way that the 
time direction is preserved. The set {ma} is a collection of non- negative integer labels attached 
to the normals a G 9ft so that raw(a) = for all w ^ Q. Huygens' principle holds for (|), (H), 
provided n is odd, and 

n > 3 + 2 ^ m„ . (7) 

In the present work we construct a new class of self-adjoint wave-type operators (||) satisfying 
Huygens' principle in Minkowski spaces. As we will see, this class provides a natural extension 
of the hierarchy of Huygens' operators associated to Coxeter groups. On the other hand, it turns 
out to be related in a surprisingly simple and fundamental way to the theory of solitons. 

To present the construction we consider a + l)-dimensional Minkowski space M""*"^ = 
Ri'" with the metric signature (+, —,—,...,—) and fixed time direction Q G M"+^. We write 
CsYxin + 1,2) C Gr(n + 1,2) for a set of all 2-dimensional space-like linear subspaces in M""^^ 
orthogonal to Q. Every 2-plane E G Gr_|_(?T, + 1, 2) is equipped with the usual Euclidean structure 
induced from M"+^. To define the potential u[x) we fix such a plane E and introduce polar 
coordinates (r, yj) therein. 

Let {ki)fLi be a strictly increasing sequence of integer positive numbersQ: < ki < k2 < . . . < 
k^^i < kjsf, and let {\E'i (</?)} be a set of 27r-periodic functions on R^: 

"ifiiip) := cos{ki(f + (fi) , G R , (8) 

associated to (ki). The Wronskian of this set 



>V[^i,^2,---,^7v] := det 



/ *i(y?) '^2{^) ■■■ ^7v(</') \ 
%{V) %{v) ... vl/^(^) 



(9) 



^Using the terminology adopted in the group representation theory we will call such integer monotonic sequences 
partitions. 
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does not vanish indentically since "^ii^f) are linearly independent. 
Let 



N 



X e M"+^ I rl'^l W [^1, ^2, . . . , ^7v] = , \k\:=J2ki 



i=l 



be an algebraic hypersurface of zeros of the Wronskian (|^) in the Minkowski space M""*"^, and let 
Q C M"^^ \ H be an open connected part in its complement. 

We define u{x) in terms of cylindrical coordinates in M""*"^ with polar components in E : 



2 / d^^ 



u = Uk{x):=--\^—j \ogW[^,{ip),^2{^),...,^N{^)] . (10) 

It is easy to see that in a standard Minkowskian coordinate chart u{x) is a real rational function 
on M"^"*"^ having its singularities on S. In particular, it is locally analytic in Vt. 
Our main result reads as follows. 

Theorem. Let M"+^ = R^'" he a Minkowski space, and let 

:= □„+! +Uk{x) (11) 



he a wave-type second order hyperholic operator with the potential (|T^ associated to an arhitrary 
strictly monotonic partition {ki) of height N : 

< ki < k2 < . . . < kN , heZ, i = l,2,...,N. 

Then operator Ci^k) satisfies Huygens' principle at every point ^ G fi, provided n is odd, and 

n>2kN + 3. (12) 



Remark I. A similar result is also valid if one takes an arbitrary Lorentzian 2-plane H G 
Gr||(n + 1,2) in the Minkowski space M"+^ containing the time-like vector 6. More precisely, in 
this case the potential Uk{x) associated to the partition (ki) is introduced in terms of pseudo-polar 
coordinates (f?, i?) in H: 

2 / \ ^ 

Uk{x) := -— ( — ] logW[V^i,^2, •••,^iv] , (13) 

where x^ = g sinh ■(9 , and, say, x^ = g cosh -(9 . The functions ipi involved in ([13|) are given by 

= cosh(fci ^ + , ^ieR . (14) 

The theorem formulated above holds when the potential (|TD[) is replaced by (|T^. 

Remark II. The potentials (|I0|) considerably extend the class of Calogero-Moser potentials 
(P) related to Coxeter groups of rank 2. Indeed, in any Coxeter group ^ is a dihedral group 
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hil), i-e. the group of symmetries of a regular 2g-polygon. It has one or two conjugacy classes 
of reflections according as q is odd or even. The corresponding potential can be rewritten in 
terms of polar coordinates as follows (see [0): 

m(m + l)g2 

sm {q(p) 

and 

m(m + l)(g/2)^ mi(mi + 1) (g/2)^ 

w(r, = — . ^ + ^ , when q even . 

r^sm {q/2)(f cos^(g/2j 

It is easy to verify that formula ([TUp boils down to these forms if we fix := m; (fi : = 
(— l)*7r/2 , i = 1, 2, . . . , , and choose 

k := (g, 2q, 3g, . . . , mq) , 

when q is odd, and 

^ '■= ^> y ' • • • ' ("^ - "^1)2' ^ + ("^ ~ "^1)2' 2g + (m - mi)-, . . . , (m + mj-) , 
when q is even and m > mi , respectively. 

Remark III. Let us set ipi = Akft + ipoi and = —Akft + 'doi , i = 1,2, . . . , N ; ipoi, 'doi E R. 
The angular parts of potentials (pUj), (P^, i.e. 

^^(¥') = -2(^) log>V[vI/i((^),vI/2(<^),...,vl>^(y,)] , (15) 



are known (see, e.g., |^) to be respectively singular periodic and proper N-soliton solutions 
of the Korteweg-de Vries equation 

vt = -v^^^ + Qvv^ . (17) 

It is also well-known that A^-soliton potentials (|16|) constitute the whole class of so-called reflec- 
tionless real potentials for the one-dimensional Schrodinger operator L = —d'^/dd'^ + v{'d) (see, 
e.g., i). 

In conclusion of this section we put forward the following conjecture]^. 



Conjecture. The wave-type operators ( pT]) with potentials of the form ([10|) give a complete 
solution of Hadamard's problem in Minkowski spaces M"+^ within a restricted class of linear 
second order hyperbolic operators 



Note added in the proof. This conjecture has been proved recently by one of the authors in 
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with real locally analytic potentials u = ^(x^jX^) depending on two spatial variables and homo- 
geneous of degree (—2): u{ax^, ax^) = a~'^u{x^, x^) , a > . 
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II. Huygens' principle and Hadamard-Riesz expansions 

The proof of the theorem stated above rests heavily on the Hadamard theory of Cauchy's problem 
for linear second order hyperbolic partial differential equations. Here, we summarize briefly some 
necessary results from this theory following essentially M. Riesz's approach |^ (see also [ |21| |, 

il). 

Let M"'^^ = R^'" be a Minkowski space, and let Q be an open connected part in M"+^. We 
consider a (formally) self-adjoint scalar wave-type operator 

£ = □„+! + m(x) , (18) 

defined in Q, the scalar field (potential) u{x) being assumed to be in C°°{Q). For any ^ G fi, we 
define a cone of isotropic (null) vectors in M"+^ with its vertex at ^: 

7(x, := (x° - er - - er - ... - (x- - cr = o , (i9) 

and single out the following sets : 

C±(0 := {x e M"+i I 7(x, = , ^0 ^ xO} , 

(20) 

J±(0 := {xGM"+M 7(a:,0 >0,e°^a;0} . 

Definition. A (forward) Riesz kernel of operator £ is a holomorphic (entire analytic) mapping 
A 1-^ $^(x,,^) , A G C, with values in the space of distributions^ T''(fi), such that for any ^ G fi: 

(i) supp$^(x,OC J+(0, 

(iz) C[<l>^{x,0]=<^Uix,0 , (21) 

(m^ ^^{x, = 5(x - . 

^By a distribution / G 'D'{il) we mean, as usual, a linear continuous form on the space 'D{n) of C°°-functions 
with supports compactly imbedded in ^l (cf., e.g., [p3|). 
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The value of the Riesz kernel ^) := ^) at A = 1 is called a (forward) fundamental 

solution of the operator C: 



£[<l>+(x, 0] = S{x - , supp <l>+(x, C J+(0 . (22) 
Such a solution is known to exist for any u{x) G C°°{Q), and it is uniquely determined. 

Definition. The operator C defined by (p!8| ) satisfies Huygens' principle in a domain fio ^ ^ 

in M"+^ if 

supp$+(x,0CC+(0=9J+(0 . (23) 

for every point ^ & Qq. 

The analytic description of singularities of Riesz kernel distributions (and, in particular, fun- 
damental solutions) for second order hyperbolic differential operators is given in terms of their 
asymptotic expansions in the vicinity of the characteristic cone by a graded scale of distribu- 
tions with weaker and weaker singularities. Such "asymptotics in smoothness", usually called 
Hadamard-Riesz expansions, turn out to be very important for testing Huygens' principle for the 
operators under consideration. 

In order to construct an appropriate scale of distributions [Riesz convolution algebra) in 
Minkowski space M"'^^ we consider (for a fixed G M""^^) a holomorphic P'- valued mapping 
C — *■ ©'(M^^^) , A I— > R\{x,^), such that R\{x,^) is an analytic continuation (in A) of the 
following (regular) distribution: 

{R,ix,0,9ix))= [ ^^"^^^''J g{x) dx , ReA>^, (24) 
J+iO 

where dx = dx^ A dx^ A ... A dx^ is a volume form in M"+^, g{x) G X'(M"'+-'^), and if„+i(A) is a 
constant given by 

/7„+i(A) = 27r V4^~ir(A)r (A - (n - l)/2) . (25) 

The following properties of this family of distributions are deduced directly from their defini- 
tion. 

For all A G C and ^ G M"+^ we have 



snppR,ix,O^J+iO (26) 

^n+i Rx = Rx-1 , (27) 

Rx*R^ = i?A+M , /i G C , (28) 

ix-^,d^)Rx = i2X-n + l)Rx , (29) 

YRx = r (A), (A - (n - l)/2)^ , z/ G Z>o , (30) 

where {k)u '■= T{k, + z/)/r(K) is Pochhammer's symbol, and 7 = 7(0;,^) is a square of the geodesic 
distance between x and in M"+^. 
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In addition, when n is odd, one can prove that 

^A(a:,0 = ^^^p^^ for A = l,2,...,(n-l)/2, (31) 



where (7) stands for the m-th derivative of Dirac's delta-measure concentrated on the surface 
of the future-directed characteristic half-cone C+(^). 

Another important property of Riesz distributions is that 

Ro{x,0 = Hx-0 ■ (32) 

Formulas (^Bf), (^71) , (|3^) show that Rx{x,C,) is a Riesz kernel for the ordinary wave operator 
□n+i. The property (^) means precisely that in even- dimensional Minkowski spaces M"'^^ {n is 
odd) Huygens' principle holds for sufficiently low powers of the wave operator , d < {n — l)/2. 

Now we are able to construct the Hadamard- Riesz expansion for the Riesz kernel of a general 
self-adjoint wave-type operator (0) on M""*"^. 

First, we have to find a sequence of two-point smooth functions Uu := Uy{x,^) G C°°{yt x 
Vt) , 1/ = 0, 1, 2 . . ., as a solution of the following transport equations: 

{x-^,d^)U,{x,0 + i'U,{^,0 = -\'^[Uu-i{x,0] , '^>1- (33) 

It is well-known (essentially due to pGf) that the differential- recurrence system (|33D has a unique 
solution provided each U,^ is required to be bounded in the vicinity of the vertex of the characteristic 
cone and Uq{x,E,) is fixed for a normalization, i.e. 

[/o(x,0 = l, f/.(e,0~(^(l) , Vz/=1,2,3,... 

These functions Ui^ are called Hadamard' s coefficients of the operator C. 

In terms of U,j the required asymptotic expansion can be presented as follows: 

00 

^"lix.O -Y.^''{>^)M^,ORx+u{x,0 . (34) 

One can prove that for a hyperbolic differential operator C with locally analytic coefficients the 
Hadamard- Riesz expansion is locally uniformly convergent. From now on we will restrict our 
consideration to this case. 

For A = 1 formula (0) provides an expansion of the fundamental solution of the operator C 
in a neighborhood of the vertex a; = of the characteristic cone: 

00 

'^'+(a;,e) = 5Z4V![/,(x,Oi?.+i(x,0 • (35) 



When n is even, we have supp Rt,j^i{x,^) = J+{^) for all = 0, 1, 2, . . ., and therefore Huygens' 
principle never occurs in odd-dimensional Minkowski spaces M^'"*"^. 
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On the other hand, in the case of an odd number of space dimensions n > 3, we know due 
to ( PTD that for 1/ = 0, 1, 2, . . . , (n — 3)/2, supp R^^i{x,^) = . Hence, using (|50D, we can 

rewrite the series (^) in following form: 

<^+ix,0 = ^^{vix,0St'\l) + Wix,0v4l)) , (36) 
where p := (n — l)/2 , r7_|_(7) is a regular distribution characteristic for the region J+(0- 

{v+i7),9ix))= j gix)dx, gix)eViM-+') , 

and V{x,C,) , W{x,C,) are analytic functions in a neighborhood of the vertex x = ^ which admit 
the following expansions therein: 

V{x, = £ ' U^x, Y , (37) 



u=0 



oo 



The function W[x^(^ is usually called a logarithmic term of the fundamental solution^. 

It follows directly from the representation formula (^) that operator L satisfies Huygens' 
principle in a neighborhood of the point ^, if and only if, the logarithmic term W[x^^ of its 
fundamental solution vanishes in this neighborhood identically in x: W[x^^ = . 

The function W{x,C,) is known to be a regular solution of the characteristic Goursat problem 
for the operator C : 

C[Wix,O] = (39) 



with a boundary value given on the cone surface C+(^) . Such a boundary problem has a unique 
solution, and hence, the necessary and sufficient condition for £ to be Huygens' operator becomes 

W{x,O=0, (40) 



where the symbol = implies that the equation in hand is satisfied only on C+(^) . By definition 
(|38|) , the latter condition is equivalent to the following one 

71 — 1 

U,{x,O = 0, P=^- (41) 

In this way, we arrive at the important criterion for the validity of Huygens' principle in terms of 
coefficients of the Hadamard-Riesz expansion (0). Equation ( PD is essentially due to Hadamard 



[26|. It will play a central role in the proof of our main theorem. 

^ Such a terminology goes back to Hadamard's book |26| , where the function W{x, ^) is introduced as a coefRcient 
under the logarithmic singularity of an elementary solution (see for details |15|, pp. 740-743). 
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III. Proof of the main theorem 

We start with some remarks concerning the properties of the one-dimensional Schrodinger operator 

/ o \ 2 

kk) ■■=-{ — ] (42) 



with a general periodic soliton potential 



dip 



Vkiif) := -2 log W [^1, . . . , ^7v] . (43) 

Here, as already discussed in the Introduction, W [^'i, ^^2, • • • , ^tv] stands for a Wronskian of the 
set of periodic functions on R^: 

"ifiiip) := cos(A;i (f + (pi) , (pi eK , (44) 

associated to an arbitrary strictly monotonic sequence of real positive numbers ("soliton ampli- 
tudes"): < fci < . . . < kN_i < kjsf . 

It is well-known (see, e.g., |^^) that any such operator Lj-^) (as well as its proper solitonic 
counterpart ([16|) ) can be constructed by a successive application of Darhoux-Crum factorization 
transformations ([]IT|, |T^) to the Schrodinger operator with the identically zero potential: 



To be precise, let L be a second order ordinary differential operator with a sufficiently smooth 
potential: 



We ask for formal factorizations of the operator 

L -\I = A* o A , (47) 

where / is an identity operator, A is a (real) constant, and A , A* are the first order operators 
adjoint to each other in a formal sense. 

According to Frobenius' theorem (see, e.g., [^), the most general factorization (|47|) is obtained 
if we take x(v^) ^-s a generic element in Ker(L — A /) \ {0} and set 

Indeed, A*oA is obviously self-adjoint second order operator with the principal part —d"^/ dip^ . 
Hence, it is of the form (^Bf). Moreover, since y4[x] = 0, we have x ^ KerA* o A , so that ( ^Tf ) 
becomes evident. 
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Note that for every A G R we actually get a one-parameter family of factorizations of L — XI. 
This follows from the fact that dim Ker(L — XI) = 2, whereas xif) cind C xif) give rise to the 
same factorization pair {A, A*) . 

By definition, the Darboux-Crum transformation maps an operator L = XI + A*oA into the 
operator 

L:=XI + AoA\ (49) 

in which A and A* are interchanged. The operator L is also a (formally) self-adjoint second-order 
differential operator 



L :-- 



(50) 



where v{(p) is given explicitly by 



v{^) = v{^) - 2 ( ^ ) log xiv) ■ 



(51) 



The initial operator L and its Darboux-Crum transform L are obviously related to each other via 
the following intertwining indentities: 



L o A = Ao L 



LoA* =A*oL . 



(52) 



The Darboux-Crum transformation has a lot of important applications in the spectral theory 
of Sturm-Liouville operators and related problems of quantum mechanics [0. In particular, it is 
used to insert or remove one eigenvalue without changing the rest of the spectrum of a Schrodinger 
operator (for details see the monograph |^ and references therein). 

The explicit construction of the family of operators (^) with periodic soliton potentials (^3[) 
is based on the following Crum's lemma: 



Lemma ( []T^ ) . Let L be a given second order Sturm-Liouville operator ( PB| ) with a sufficiently 
smooth potential, and let {\E'i, ^^2, • • • , ^Af} be its eigenfunctions corresponding to arbitrarily 
fixed pairwise different eigenvalues {Ai, A2, . . . , Xn} , i.e. e Ker(L — Xil) , i = 1,2, . . . , N . 
Then, for arbitrary \1' G Ker(L — A /) , A G R , the function 



>V[^l,^2,---,^jV,^] 
W[^'l,^'2,...,^'iv] 



(53) 



satisfies the differential equation 



with the potential 



d 



Xiv(<^) = ^Xn{v) 



Vn{^) := v{ip) - 2 ( ^ ) log W [^1, ^^2, • • • , "^Nl 



(54) 



(55) 
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Given a sequence of real positive numbers {ki)^i: < ki < k2 < ■ ■ ■ < , the Darboux-Crum 
factorization scheme: 



Li := Ai.i o A*_^ + kjl = A* oAi + kj^^ I L^+i := AiO A* 
starting from the Schrodinger operator ( ^Sj ) with a zero potential 



k'^ T 



(56) 



A) =AloAo + kll 



produces the required operator L(fc) = with the general periodic potential (^). 

Now we proceed to the proof of our main theorem formulated in the Introduction. 

When = , the statement of the theorem is evident, since the operator is just the 
ordinary wave operator in an odd number n of spatial variables. 

Using the Darboux-Crum scheme as outlined above we will carry out the proof by induction 
in N. 

Suppose that the statement of the theorem is valid for all m = 0,1,2, N. Consider an 
arbitrary integer monotonic partition (ki) of height : < ki < k2 < ■ ■ ■ < kN , G Z. 
By our assumption, the wave-type operator 



N 



(fc) 



□„+i + Uk{x) 



(57) 



associated to this partition, satisfies Huygens' principle in the (?2+l)-dimensional Minkowski space 
M"^^ with n odd, and n > 2 fc^v + 3 . We fix the minimal admissible number of space variables, 
i.e. n = 2 /ctv + 3 , and denote 

7X — 1 

P ■■= = kN + l. (58) 

By construction, the operator Cn can be written explicitly in terms of suitably chosen cylin- 
drical coordinates in M"+^: 



C 



N 



□ 



n-1 



9_ 

dr 



+ 



ld_ 

r dr 



d_ 

dip 



(59) 



where (r, (p) are the polar coordinates in some Euclidean 2-plane E orthogonal to the time direction 
in M'^^^, i.e. E G Gr_L(n + 1,2); is a wave operator in the orthogonal complement 

E^ = M"-i of E in M"+i; and vnI'p) is a 27r-periodic potential given by (H). 
Let k := kj^^i be an arbitrary positive integer such that 



k > k 



N 



(60) 



We apply the Darboux-Crum transformation ( |56D with the spectral parameter k to the angular 
part of the Laplacian in E. For this we rewrite in the form 



■■N 



□n-1 - 



(61) 
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and set 



dV ld__l_ 
dr j r dr 



{An oA% + k^) 



(62) 



where := An^^p) and A*^ := A*^{({)) are the first order ordinary differential operators of the 
form (gSp. 

According to (p2D, we have 

Cn+1 o An = An o Cn , Cn o A*^ = A*^ o Cn+i ■ (63) 

Let $^(x,^) and $^^^(x,^) be the Riesz kernels of hyperbohc operators Cn and Cn+i 
respectively. Then, by virtue of (p3| ) we must have the relation 

AUv) [$f +1] - An{<P) [$f ] =0 for all A G C , (64) 

where AAr(0) is the differential operator An written in terms of the variable conjugated to 
. Indeed, if identity (|6^) were not valid, one could define a holomorphic mapping $^ : C — > 
P', A (a;,0 , such that 

^^(^,0 :=$^(x,0 + a ^(y') - A;v(0) . (65) 

The distribution $^(a;,^), depending on an arbitrary complex parameter a G C, would also 
satisfy all the axioms (|2l| ) in the definition of a Riesz kernel for the operator Cn- In this way, we 
would arrive at the contradiction with the uniqueness of such a kernel. 

In particular, when A = 1, the identity (|^) gives the relation between the fundamental solu- 
tions $^(x, = <l>f(x,0 and <l>++^(a;,0 = <l>f+^(x,0 of operators Cn and Cn+i- In accor- 
dance with (|36D , we have 

^+{^^0 = ^ {vNix,05t'\l) + WN{x,0v+il)) (66) 

and 

<^>r'(^,0 = ^ (v^+i(^,0^r'^(7) + W^^+i(a:,0^+(7)) , (67) 

where 7 is a square of the geodesic distance between the points x and ^ in M"+^. Substituting 
(|66|), ( |67D into (|6lD, we get the relation between the logarithmic terms Wn{x,^) and Wn+i{x,^) 
of operators and Cn+i 

A%{^) [Wn+i{x,0] - An{^) [Wn{x,0] = . (68) 

By our assumption, Cn is a Huygens' operator in M""*"^, so that Wn{x,^) = 0. Hence, equation 
(|68|) implies A*^{{p) [Wn+i{x,^)] = . On the other hand, as discussed in Sect. II, the logarithmic 
term Wn+i{x,C,) is a regular solution of the characteristic Goursat problem for Cn+i, i-e. in 
particular, 

Cn+1 [Wn+i{x,O] = - (69) 
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Taking into account definition (^) of tlie operator Cn+i , we arrive at tlie following equation for 



According to (|38|) , the logarithmic term PVjv+i admits the following expansion 



(70) 



7 



v—p 



P 



n — 1 



(71) 



where U^{x,^) are the Hadamard coefficients of the operator Cn+i ■ Since the potential of the 
wave-type operator ^at+i depends only on the variables r, ip , its Hadamard coefficients Ui, must 
depend on the same variables r, ip and their conjugates p, (f) only: 



for all z/ = 0,1,2, 



(72) 



This follows immediately from the uniqueness of solution of Hadamard's transport equations (p3D. 
On the other hand, since 



7 = — — + 2r p cos{(p — (j)) 



(73) 



where s is a geodesic distance in the space E-^ = M"^^ orthogonally complementary to the 2-plane 
E, we conclude that Wn+i is actually a function of five variables: Wn+i = PVAr+i(s, r, p, 0) . 
On the space of such functions the wave operator in E-^ acts in the same way as its "radial 

part", i.e. 



— Wn+1 

S OS 



Hence, equation (|70D becomes 

dr 



d_ 

ds 



n-2 d Id k\ 

TT + - 7^ 2 I ^^+1 = • 

s OS r or 



(74) 



Now we substitute the expansion (|71|) 

oo 

Wn+1 = ^ U^{r,(p,p,< 



1 



v—p 



n — 1 



v=p 



P 



(75) 



into the left-hand side of the latter equation and develop the result into the similar power series 



in 7, taking into account formula (ffSl). After simple calculations we obtain 



E 



v=p 



1 k"^ 



4 (r - p cos(y? - (j))) U'^^^- 



(76) 
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- 2 (^2 (z/ + 1) - ^ cos(<^ - 0)) f/,+1 - 4p2 sin^if - 0) 



7 



u—p 







where the prime means differentiation with respect to r . 

Since the functions f/^ do not depend exphcitly on 7 , equation (|76|) can be satisfied only if 
each coefficient under the powers of 7 vanishes separately. In this way we arrive at the following 
differential-recurrence relation for the Hadamard coefficients of the operator C^^i : 



4p2 sin2(^-0)f/^+2 



U" + -U' 
r 



2p 



+ cos((^ - 0) (2 r f/;+i + f/,+i) - 4 (r f/^+, + + U,+i) 



(77) 



where runs from p: u = p, p + 1, p + 2, . . . 

To get a further simplification of equation ([77| ) we notice that all the Hadamard coefficients 



of the operators under consideration ([Tl|) , (|ToD are homogeneous functions of appropriate degrees. 
More precisely, they have the following specific form 



(rpy 



(T^(y?,0) , i/ = 0,l,2,... 



(78) 



where (Tu{<f, <p) = o'ui4>, f) are symmetric 27r-periodic functions depending on the angular variables 
only. 

In order to prove Ansatz ([78|) we have to go back to the relation (|6^) between the Riesz kernels 
of operators Cn and ^at+i : 



A*M [<|.f +i(x, 0] - A^{4>) [$f (x, 0] = , A G C 



(79) 



If we substitute the Hadamard- Riesz expansions (0) of the kernels $^(x, ^) and $j^'*'^(x, ^) 
into ([79|) directly and take into account that and its adjoint A*j^ are the first order ordinary 
differential operators of the following form (cf. (^H])): 



d d 

where fN^f) = {9/ dip) logXAr(v^) , we obtain 

00 

^ 4^^ (A). \2rp sin(^ - 0) {u!!,f - f/.^+J - 



^0) 



u=0 



d 



d 



R 



where {r,ip, p, cj)) and U^~^^{r,{p, p,(j)) are the Hadamard coefficients of operators Cn and 
Cn+1 respectively; R\ := R\{x,^) is the family of Riesz distributions in M"+^ . 
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The same argument as above (see the remark before formula ([77D) shows that all the coefficients 
of the series (^TJ) under the Riesz distributions of different weights must vanish separately. So we 
arrive at the recurrence relation between the sequences of Hadamard's coefficients of operators 
C]^ and >Cjv+i • 



TjN+l _ TjN 1 



2rp sm{ip — 0) 



u 



(82) 



where U^^'^ = Uq = 1 and = 0, 1, 2, . . . Now it is easy to conclude from (^) by induction in 
that the Ansatz ( ffSD really holds for Hadamard's coefficients of all wave-type operators (|l5) with 
potentials ([T0|) . 

Returning to equation (|77| ) and substituting (|78D therein, we obtain the following three-term 
recurrence relation for the angular functions cri,{ip, (p) : 

4 sm^{ip - (p) au+2 = {v"^ - k^) (Ju - 2(2z/ + 1) cos(v9 - ^) a^+i , (83) 

where u = p, p + 1, p + 2, . . . . 

In order to analyze equation (|83D it is convenient to introduce a formal generating function for 
the quantities {cu} '■ 

m--=Y,^A^A)j^,- (84) 

The recurrence relation (^3[) turns out to be equivalent to the classical hypergeometric differ- 
ential equation for the function F{t) 

d'^F rIF 
(4(l-cu2)+4cut-t2) __ + (2p+l)(2cu-t) — + (k^-p^)F = 0, (85) 



where u := cos{(p — 0) . The general solution to ( ^5]) is given in terms of Gauss' hypergeometric 
series: 

F{t) = C 2F^{p - k;p + k-p+ 1/2 \z) + Ci z'^^'^^ 2 Fi(l/2 - k- 1/2 + k-3/2-p\z) , (86) 
where z := (t — 2c<j + 2)/4 and 2 Fi is defined by 

2F,{a;b-c\z):=}_^ f ^ - . (87) 

As discussed in Sect. II, the Hadamard coefficients U^{x,^) must be regular in a neighborhood 
of the vertex of the characteristic cone x = C, . When x ^ ,^ , we have uj ^ 1 and Up{^,^) oc 
crp(0, 0) = F{0)\i^=i is not bounded unless Ci = . 

In this way, setting Ci = in (|86D , we obtain 



£ ^-('/'^ ^) jfl-y^ =C 2F,ip-k;p + k;p+l/2\it-2uj + 2)/4) . 
i/=p ^ ^' ' 



16 



Now it remains to recall that by our assumption (60) k G Z and k > . Since p = (n — l)/2 
kj\f + 1 , we have k > p . So the hypergeometric series in the right-hand side of equation 
is truncated. In fact, the generating function (|84D is expressed in terms of the classical Jacobi 
polynomial 'P^^Zp^'^'^^^^'^\uJ — t/2) of degree k — p. Hence, ak+i{f,(f)) = , and the {k + l)-th 
Hadamard coefficient of the operator Cn+i vanishes identically: 

Uk+i{x,O = 0. (89) 

According to Hadamard's criterion (|41|), it means that the operator Cn+i satisfies Huygens' prin- 
ciple in Minkowski space M"+^ , if n is odd and 

n>2k+3 , 

Thus, the proof of the theorem is completed. 



IV. Concluding remarks and examples 

In the present paper we have constructed a new hierarchy of Huygens' operators in higher di- 
mensional Minkowski spaces NP~^^,n > 3. However, the problem of complete description of the 
whole class of such operators for arbitrary n still remains open. As mentioned in the Introduction, 
the famous Hadamard's conjecture claiming that any Huygens' operator C can be reduced to 
the ordinary d'Alembertian with the help of trivial transformations is valid only in M'^"''^ . 

Recently, in the work |Q one of the authors put forward the relevant modification of Hadamard's 
conjecture for Minkowski spaces of arbitrary dimensions. Here we recall and discuss briefiy this 
statement. 

Let Q be an open set in Minkowski space M"+^ = R"'"'"^ , and let J-'{Q) be a ring of partial 
differential operators defined over the function space C°°{Q) . For a fixed pair of operators Cq, C G 
jF(fi) we introduce the map 

^dc,Co ■■ H^) ^ H^) , A ^ ^dc,Co[A] , (90) 

such that 

adc,Co[A] ■.= CoA-AoCo. (91) 
Then, given M G Z>o , the iterated ad£,£o-map is determined by 

^dZco [A] := e.dc,Co [ad^,£o [• • • ^dc,Co [A]] . . ] = J^i-l)' , oAoC',. (92) 

k=o ^ ^ 

Definition. The operator C G J-'{Q) is called M-gauge related to the operator Co G !F{fl) , if 
there exists a smooth function 6{x) G C°^{Q) non- vanishing in Q, and an integer positive number 
M G Z>o , such that 

^dc,Co [^(^)] = identically in J^{n) . (93) 
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In particular, when M = 1 , the operators L and £o are connected just by the trivial gauge 
transformation L = 9{x) o Cq o 6{x)^^ . 

The modified Hadamard's conjecture claims: 

Any Huygens' operator C of the general form 

C = n^+i + {a{x),d)+u{x) , (94) 

in a Minkowski space M""*"^ (n is odd, n > 3) is M -gauge related to the ordinary wave 
operator □„+! in M*^^^. 

For Huygens' operators associated to the rational solutions of the KdV-equation (^, (|^) and to 
Coxeter groups (|), (|]) this conjecture has been proved in and 0. In these cases the required 
identities (0) are the following 



<^:[n(^°)]=0, M,:=M^, (95) 
where is given by (H) with the potential (H) for /c = 0, 1, 2, . . . and 

adfr,SK(^)] = ' := ^ m« , (96) 

where Cm is defined by (|), (g) and 7rm(x) := naG5R+ ^)"'" • 

It is remarkable that for the operators constructed in the present work the modified Hadamard's 
conjecture is also verified. More precisely, for a given wave-type operator 

2 f d^^ 



Cik) = Dn+l - ^ i^g^j logW [^l(^), ^2(<^), . . . , ^Niv)] , (97) 

associated to a positive integer partition (hi) : < ki < k2 < ■ ■ ■ < , we have the identity 

^4l!'co^Qik){x)]=0, (98) 

where 9(fc)(x) := r''^' W ^^2, • • • , ^tv] and \k\ := ki is a weight of the partition (ki) . 

We are not going to prove (|98D in the present paper. A more detailed discussion of this 
identity and associated algebraic structures will be the subject of our subsequent work. Here, we 
only mention that such type identities naturally appear [ffl-p[ in connection with a classification 



of overcomplete commutative rings of partial differential operators |13|, [Q, [41|, and with the 



bispectral problem |20 



We conclude the paper with several concrete examples illustrating our main theorem. 
1. As a first example we consider the dihedral group hiq), Q G Z>o, acting on the Euclidean 
plane E = C Gr_|_(n + 1,2) and fix the simplest partition k = (q) and the phase ip = tc/2. 
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According to Remark II, in this case our theorem gives the wave-type operator with the Calogero- 
Moser potential related to the Coxeter group hiq) with m = 1: 

r^sm {q(p) 

This operator satisfies Huygens' principle in M"+^ if n is odd and n > 2 g + 3. The Hadamard 
coefficients of C(^k) can be presented in a simple closed form in terms of polar coordinates on E: 

Uo = l, 

^r^ 1 Tj-HcOs(v. - 0)) ^ ^^^^ 

(2rpY sin(g(y9) sin(g0) ' ~ ' 

where Tq{z) := cos(g arccos(2;)), z G [—1,1], is the g-th Chebyshev polynomial, and Tq'^\z) is 
its derivative of order u with respect to z. These formulas are easily obtained with the help of 
recurrence relation (^21). 

2. Now we fix iV = 2, ki = 2, k2 = 3 and ipi = n/2, ip2 = ■ The corresponding wave-type 
operator 

, 10 {xl + xl){15xl-xl) 

(5X2 + ^l) ^1 

satisfies Huygens' principle for odd n > 9. The nonzero Hadamard coefficients of this operator 
are given explicitly by the formulas: 

?7o = l, 

40a;2ei6a;i + l^iix^" + 75 ^'^a' + 15 ^'a:!^ -h^^^Xi^ 



U2-- 



2eiXi(5X22+Xi2)(5e2' + ei') 

120x2^162:1 + 15 6'2:i2 -bii^xi^ + lbii^X2^ + 75 6=^X22 



46V(5x22 + xi2) (56' + 6') 
15 X26 

~eiVM5x22 + xi2) (5 6' + 6') 



3. Now we take = 3, the partition k = (1, 3, 4) , and the phases = tf2 = = 7i"/2 . The 
corresponding operator 

_ , 12 (49xf + 28x1x1 - 4) 

J~-{k) - ^n+l H 2 2 , 2^2 ' 

x^ (7xf + X^) 

is a Huygens operator in M"+^ when n is odd and n > 11 . The nonzero Hadamard's coefficients 
are 

Uo = l, 

-21 ^'xi^ - 42x2ei6a;i - 2Ui'a;2' + 3^2^X2'' - U7 ^i^i^ 



6a:2(7xi2 + x22) (76' + 6') 
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_ 735 ^I'xi^ + 504 S266a;i + 105 ^1^2' - 21 ^'xa^ + 105 ^'xi^ 

(72 — 



4e2V(7a:i2 + x22)(7ei' + 6') 
-1260 a;2^i62;i + 21 ^2^X2" - 105 ^'2:2' - 105 ^'^i^ - 735 ^^i^ 



86V(7xi2 + x22) (76' + 6') 
315 



46V(7xi2 + x22)(7ei' + 6') ■ 

4. The last example illustrates Remark I following the theorem (see Introduction). In this 
case we consider the operator (|TT]) with the potential (|TB|) associated with the proper A^-soliton 
solution of the KdV equation. We take N = 2 and fix A;i = 1, A;2 = 2 . The real phases are chosen 
as follows = arctanh (1/2), i)2 = arctanh(l/4) . The corresponding operator C(^k) reads 

2 (2x0 - 3x1) {3x1 - Q^oxj + 4:Xixl + 8x1) 



x\ (4xo — 2xqXi — x'^'^ 



According to the theorem, it is huygensian provided n is odd and n > 7. The nonzero Hadamard 
coefficients are given by the following formulas: 



Ui 



Uo = l, 

ieoxf + 9e?x2 - 16^0^x2 + 8^0^x0x1 



U2 



2x1 (4x2 - 2xoXi - X?) 6 (4^0^ - 2^06 - e?) 

8^063^0 - mfxQXi + A^iXl - U^O^ixl + 16^063^03^1 

2x1 (4x2 - 2xoXi - X?) 6 (46' - 266 - 6') 

^ 5 (2 6-6)(2xo-xi) 

4 xi (4 xo2 - 2 xoxi - xi2) 6 (4 6' - 2 66 - 6') 



References 

[1] M. J. Ablowitz, and P. A. Clarkson, Solitons, Nonlinear Evolution Equations and Inverse 
Scattering, Cambridge University Press, Cambridge, 1991. 

[2] M. Adler, and J. Moser, On a class of polynomials connected with the Korteweg-de Vries 
Equation, Comm. Math. Phys., 61 (1978), 1-30. 

[3] H. Airault, H. McKean, and J. Moser, Rational and elliptic solutions of the KdV-equation 
and related many body problems, Comm. Pure Appl. Math., 30 (1977), 95-148. 

[4] Yu. Yu. Berest, Hierarchies of Huy gens' operators and Hadamard's conjecture, preprint CRM- 
2296, to be published in Acta Math Appl. (1997). 



20 



[5] Yu. Yu. Berest, Lacunae of hyperbolic Riesz kernels and commutative rings of partial differ- 
ential operators, preprint CRM-2325 , to be published in Lett. Math. Phys. (1997). 

[6] Yu. Yu. Berest, Solution of a restricted Hadamard's problem in Minkowski spaces, preprint 
CRM-2359, to be pubhshed in Comm. Pure Appl. Math. (1997). 

[7] Yu. Yu. Berest, and Yu. A. Molchanov, Fundamental solutions for partial differential equations 
with reflection group invariance, J. Math. Phys., 36(8) (1995), 4324-4339. 

[8] Yu. Yu. Berest, and A. P. Veselov, Huygens' principle and integrability, Uspekhi Mat. Nauk, 
49(6) (1994), 5-77. 

[9] Yu. Yu. Berest. and A. P. Veselov, Huygens' principle and Coxeter groups, Uspekhi Mat. 
Nauk, 48(2) (1993), 181-182. 

[10] Yu. Yu. Berest, and A. P. Veselov, Hadamard's problem and Coxeter groups: new examples 
of Huygens' equations, Funct. Anal. Appl., 28(1) (1994), 3-15. 

[11] J.L.Burchnall, and T.W.Chaundy, A set of differential equations which can be solved by poly- 
nomials, Proc. London Math. Soc., 30 (1929-30), 401-414. 

[12] J. Carminati, and R. G. McLenaghan, On the validity of Huygens' principle for the wave 
equation on curved space-time in "Huygens' principle 1690-1990: Theory and Application^'' 
(H.Blok, H.A.Ferwerda, H.K.Kuiken, Eds.), 435-442, Elsevier Publishers B.V.,1992. 

[13] O. A. Chalykh, and A. P. Veselov, Commutative rings of partial differential operators and 
Lie algebras. Comm. Math. Phys., 126 (1990), 597-611. 

[14] O. A. Chalykh, and A. P. Veselov, Integrability in the theory of the Schrddinger operator and 
harmonic analysis, Comm. Math. Phys., 152 (1993), 29-40. 

[15] R. Courant, and D. Hilbert, Methods of Mathematical Physics U, Interscience Publ., N.-Y., 
1962. 

[16] M. Crum, Associated Sturm- Liouville systems. Quart. J. Math., 6(2) (1955), 121-127. 

[17] G. Darboux, Sur la representation spherique des surfaces, Compt. Rend., 94 (1882), 1343- 
1345. 

[18] P. A. Deift, Application of a commutation formula, Duke Math. Journal, 45(2) (1978), 267- 
310. 

[19] J. J. Duistermaat, Huygens' principle for linear partial differential equations in "Huygens' 
principle 1690-1990: Theory and Applications^' (H.Blok, H.A.Ferwerda, H.K.Kuiken, Eds.), 
pp. 273-297, Elsevier Publishers B.V., 1992. 



21 



[20] J. J. Duistermaat, and F. A. Griinbaum, Differential equations in the spectral parameter, 
Comm. Math. Phys., 103 (1986), 177-240. 

[21] F. G. Friedlander, The Wave Equation on a Curved Space-time, Cambridge University Press, 
Cambridge, 1975. 

[22] A. M. Gabrielov, and V. P. Palamodov, Huygens' principle and its generalizations in 
"/. G. Petrovskii, Collected papers" , v.l, 449-456, Nauka, Moscow, 1986 (in Russian). 

[23] I. M. Gel'fand, and G. E. Shilov, Generalized Functions, V.l, Acad. Press, N.-Y., 1964. 

[24] P. Giinther, Huygens' Principle and Hyperbolic Equations, Acad. Press, Boston, 1988. 

[25] P. Giintlicr, Huygens' principle and Hadamard's conjecture. The Mathematical Intelligencer, 
13(2) (1991), 56-63. 

[26] J. Hadamard, Lectures on Cauchy's Problem in Linear Partial Differential Equations, Yale 
University Press, New Haven, 1923. 

[27] S. Hclgason, Wave equations on homogeneous spaces, Lect. Notes in Math. 1077, Springer- 
Verlag, Berlin, 252-287 (1984). 

[28] N. H. Ibragimov, Transformation Groups Applied to Mathematical Physics, Reidel, Dordrecht, 
1985. 

[29] E. L. Ince, Ordinary Differential Equations, Dover, New- York, 1956. 

[30] L.Infeld, and T.E.Hull, The factorization method. Rev. Mod. Phys., 23(1) (1951), 21-70. 

[31] J. E. Lagnese, and K. L. Stcllmacher, A method of generating classes of Huygens' operators, 
J. Math. Mech., 17(5) (1967), 461-472. 

[32] J. E. Lagnese, A solution of Hadamard's problem for a restricted class of operators, Proc. 
Amer. Math. Soc, 19 (1968), 981-988. 

[33] M. Mathisson, Le probleme de M. Hadamard relatif a la diffusion des ondes. Acta Math., 71 
(1939), 249-282. 

[34] V. B. Matveev, and M. A. Salle, Darboux Transformations and Solitons, Springer- Verlag, 
Berlin-Heidelberg, 1991. 

[35] R. G. McLenaghan, Huygens' principle, Ann. Inst. Henri Poincare, 37(3) (1982), 211-236. 

[36] M. A. Olshanetsky, and A. M. Perelomov, Quantum integrable systems related to Lie algebras, 
Phys. Reps. 94 (1983), 313-404. 

[37] M. Riesz, L'integrale de Riemann-Liouville et le probleme de Cauchy, Acta Math. 81 (1949), 
1-223. 



22 



[38] R. Schimming, Korteweg-de Vries-Hierarchie und Huygenssches Prinzip, Dresdener Seminar 
fiir Theoretische Physik, Sitzungsberichte, N 26 (1986). 

[39] R. Schimming, Laplace-like linear differential operators with a logarithmic-free elementary 
solution, Math. Nachr., 148 (1990), 145-174. 

[40] K. Stellmacher, Ein Beispiel einer Huygensschen Diff'erentialgleichung, Nachr. Akad. Wiss. 
Gottingen, Math.-Phys. Kl, Ila, Bd. 10 (1953), 133-138. 

[41] A. P. Veselov, K. L. Styrkas, and O. A. Chalykh, Algebraic integrability for the Schrddinger 
operator and reflection groups, Theor. Math. Phys., 94 (1993), 253-275. 

[42] A. P. Veselov, M. V. Feigin, and O. A. Chalykh. New integrable deformations of quantum 
Calogero-Moser problem, Usp. Mat. Nauk., 51(3) (1996), 185-186. 



23 



